High-Tc-Superconductivity and Shadow State Formation 
in YBa2Cu306+5 and Bi2Sr2CaCu208+5 
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The normal and superconducting state of YBa2Cu306+« and Bi2Sr2CaCu208+« are investigated 
by using the mono- and bilayer Hubbard model within the fluctuation exchange approximation 
and a proper description of the Fermi surface topology. The inter- and intra-layer interactions, 
the renormahzation of the bilayer splitting and the formation of shadow bands are investigated in 
detaih Although the shadow states are not visible in the monolayer, we find that the additional 
correlations in bilayers boost the shadow state intensity and will lead to their observabihty. In 
the superconducting state we find a d^-2_yi symmetry of the order parameter and demonstrate the 
importance of inter-plane Copper pairing. 



Central questions in the theory of the high-Tc super- 
conductors are related to the importance of the multi- 
ple Cu02 layers within the High-Tc superconductors like 
BisSrsCaCusOg+a (BSCCO) or YBaaCuaOe-na (YBCO) 
and its influence on the superconducting pairing symme- 
try. In addition, the detailed shape of the Fermi surface 
(FS) is believed to be of importance for a quantitative 
description of transport and photoemission experiments 
and also for the material dependence of the transition 
temperature Tc Qj^- For bilayer cuprates, neutron scat- 
tering experiments found indications that there is an an- 
tiferromagnetic coupling between nearest-neighbor lay- 
ers in YBCO l^,^- Furthermore, recent angular resolved 
photoemission (ARPES) experiments found evidence for 
two separated bands in YBCO (|J^, that might be re- 
lated to the existence of two CuO-z bands caused by a 
inter-plane quasi particle transfer. However, the small 
experimentally observed bilayer splitting in YBCO and 
the difficulty to resolve two Cu02 bands in BSCCO §1] 
support the idea that the strong short ranged antifer- 
romagnetic order in the cuprates reduces the inter-layer 
hopping and might be responsible for the observation of 
shadows of the FS §,|o| in BSCCO. 

In this paper we study the bilayer Hubbard Hamilto- 
nian within the fluctuation exchange (FLEX) approxima- 
tion ||l^ and focus our attention in particular on YBCO 
and BSCCO compounds. We find that the antiferromag- 
netically correlated planes yield strong deformations of 
the quasi particle dispersions and suppression of the bi- 
layer splitting. In addition we observe that the shadows 
of the FS occur only when the inter-plane coupling is con- 
sidered. The superconducting state is investigated within 
the framework of the Eliashberg theory where the order 
parameter is found to have a dj.2_y2 symmetry with inter- 
and intra-layer Cooper pair formation. 

Our theory is based on the general multilayer Hubbard 
model that will be later on specified for a bilayer system: 



where the hopping integrals t 



persion ej^ in 2D k-space, 

site (layer) indices, (5*';, the Kronecker symbol, U 
local Coulomb repulsion and fj, the chemical potential. 
For the bilayer the interaction-free contribution of the 
Hamiltonian can be diagonalized with the unitary trans- 
formation U yielding an antibonding (— ) and a bonding 
band (-I-) with bare dispersion e^. Assuming that this 
symmetry holds also in the full interacting case, one can 
define for the normal state corresponding Greens func- 
tion G±(k, zwrri) with fermionic Matsubara frequencies 
LUm = (2?n + l)7rr and temperature T. 

In the superconducting state it is helpful to investi- 
gate the allowed symmetries of the order parameter in a 
multilayer system. Here we consider only singlet super- 
conductivity and for simplicity the bilayers case although 
similar symmetries apply to the multilayer Hamiltonian. 
Note that this discussion is only related to the inter-layer 
effects and not restricted to a certain in plane symme- 
try {d- or s-wave). At first by interchanging the two 
electrons of the Cooper pair, it follows for the gap func- 
tion in the layer representation that A;p(k) = A/'/(— k) 
due to the Pauli principle. Note that the frequency in- 
dices have been omitted for clarity. In addition to the in- 
plane symmetries, present for a single layer, the Hamilto- 
nian is invariant with respect to the inversion symmetry 
with inversion center between the layers: therefore, the 
gap function has a given parity P — ±1 and it follows 
Aii(k) = PA22(-k) and Ai2(k) = PA2i(-k). Using 
the in-plane symmetry k ^ — k it follows from this con- 
siderations that the gap-function for even parity pairing 
P = -|-1 is given in the layer representation by: 
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Here, the gap-function can be diagonalized by the trans- 
formation lA of the normal state and intra-band pairing 
(simultaneous intra- and interlayer pairing) occurs. For 
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odd parity pairing P 
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and only intra-layer pairing occurs. Hence the gap func- 
tion is off diagonal in the eigenvalue representation lead- 
ing to an interband pairing state. Consequently, the cor- 
responding Eliashberg equations decouple and the sym- 
metry with the larger Tc will determine the supercon- 
ducting state. However one might argue that there will 
be a change from intra- to inter-band pairing or vise versa 
for T < Tc which would be related to a second phase 
transition below Tc- So far this has not been observed 
experimentally, and we will only consider the solution 
with the larger energy gain of the condensate. Now, by 
looking at the inter-band pairing, for a given momenta 
two electrons from different bands that form a Cooper 
pair cannot origin both from the FS due the structure 
of the bonding and antibonding band and due to the bi- 
layer splitting. Hence, we expect intra-band pairing to 
be most dominant, e.g. yielding the largest energy gain 
of the superconducting phase which will be treated in the 
following. Note that this point is also supported by weak 
coupling approaches as shown by Maly et al. 

Now to treat the superconducting state of the 
multiband Hubbard Hamiltonian we use the Nambu- 
Eliashberg approach where the Green's function and the 
off-diagonal Green's function that describes supercon- 
ducting correlations are written in terms of 2 x 2 matri- 
ces. Here we have to account for an additional index for 
the layer or the eigenvalue representation . Conse- 

quently the Eliashberg equations for intra-band pairing 
can be written as follows by expanding the matrix self 
energy in terms of Pauli matrices: 



*k,A(«t^n) = T7 V ^ 7^ r ^ 

^k.A(«tJ„ = T7 V j:. T. ^ 
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x(£k'A' + Xk',A'(«^n')) 7 
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where 

-Dk',A'(«^n') = («t^n'^k',A'(j^n'))^ 



Here the term U5\^x' accounts for the Hartree contribu- 
tion which is for the diagonal elements absorbed in the 
chemical potential and £^ is the free dispersion that de- 
termines the FS shape. The expansion coefRcients of the 
diagonal self energy are iujn{^ — Z]^,x{'i^n)) and Xk.xii^n) , 
whereas 0k,A(*'^n) — ^)i,\{i^n)Z]ii.\{'i^n) is the coeffi- 
cient of the off-diagonal self energy which signals super- 
conductivity and Ak,A(«w„) is the gap function. The in- 
teractions V^'"^ and V^^'^ for the bilayer systems can 
be obtained by performing the summation of the FLEX 
diagrams |TT| in the layer representation and are simi- 
lar to the monolayer case when one considers the fact 
that a third dimension with two momenta points tt and 
is introduced. Furthermore for the bilayer system they 
have an inter-plane, V-^{iWm)i and an in-plane, V^' [ium), 
contribution: 



V++{]<i,lLO^) - 1/2 {VliiuJn) + V^{i0Jn)) 
V+'{\i,lLO,n) - 1/2 {VliiUn) - V^{iWn)) 
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-(£k',A' - v'^k',A'(«W„'))^ - $k',A'(*^n')^- 



This set of coupled Eliashberg equations is solved self- 
consistently on the real frequency axis ]l5[ |. Since most 
photoemission experiments were performed on YBCO 
and in particular BSCCO systems, it is necessary for a 
detailed comparison of our theory with experiments to 
use dispersions and Fermi surfaces that are closed around 
the (tt, tt) point. Thus the FS topology is characterized 
by 

= — [2t(cos(A:2;) + cos{ky)) + At' cos{kx) cos{ky) 
+2t" {cos{2kx) + cos(2fcy)) ± 

with intra-plane hopping integrals t = 0.25 eV, t' — 
-0.38t, t" = -0.06t and t±_ = 0.4< ^ 100 meV (l|] 
as explicit model for YBCO. Due to the similarity of 
the YBCO and BSCCO FS, we used for the later one 
the bare intra-plane dispersion of YBCO and i±(k) — 
1/4 i_L (cos(A:2;) -cos(fcy))^ with t±_ — OAt. The resulting 
FS is very similar to the experiments ||^ and to the band 
structure calculations For comparison with previ- 

ous results we take U = At but notice that we find no 
significant changes in our data up to values of J7 = 6t. 

In Fig. 1 we present out data for the effective interac- 
tion Vq^Lu) in the two dimensional Brillouin zone (BZ) for 
a; = where we have for simplicity neglected the bilayer 
coupling {t± — 0). Note that a; = 1 — n is the dop- 
ing concentration whereby n is the occupation number 
per site. Interestingly, we find commensurate peaks in 
the effective interaction for this compound as observed 
in neutron scattering experiments This results is 

caused by the electronic correlations since for J7 = 
pronounced incommensurabilities are present due to the 
fact that 2 kps 7^ (tt, tt) with Fermi surface momentum 
kps that shift the peaks in Vq{uj) away from Q. How- 
ever for finite U and strong scattering rates these incom- 
mensurable structures are smeared out due to the over- 
damped nature of the spin excitations and a peak at Q 
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remains. So far only the in-plane correlations in an in- 
sulated CM02-plane are considered. However, when var- 
ious layers are coupled like in the bilayer systems YBCO 
and BSCCO the influence of the inter-plane interaction 
V^^{lj) on the superconducting state and Tc is a sig- 
nificant and important open question. By investigating 
Vq^ {lu) for the YBCO and BSCCO parameterization with 

t± = 0.4t we find that the commcnsurabillitics in Vq {uj) 
remain or get even more pronounced while V^{uj) has 

the same k-dependence as Vq (w) but with an opposite 
sign. This yields an antiferromagnetic correlated bilayer 
via the hopping t± that is strongest for low doping con- 
centrations, where inter- and intra-layer interactions be- 
come comparable in magnitude although we always find 
V^^{u;)<V^{c.). 

The formation of shadow states is an important issue 
in the recent research, since it concerns the transition of 
the quasi particle excitations with doping from the sim- 
ple paramagnetic state for large x to the antiferromagnet 
for low X and from large FS for optimally doped systems 
to small FS hole pockets for low doping. Recently we 
discussed for a simple model dispersion with t' = t" = 
that the dynamical antiferromagnetic short range order 
in the cuprates leads to a transfer of spectral weight from 
the FS at kps to its shadow at kps + Q |lq | that might 
lead to the observed ARPES spectra [p|,[lO[. In this con- 
text the importance of the quasi two dimensional nature 
of the magnetic excitations for shadow states has been 
discussed recently in Ref . . Now, in Fig. 2 we demon- 
strate the influence of the realistic YBCO or BSCCO dis- 
persion on the shadow states for t± = 0. Here we plot the 
spectral density Pk(^) for two doping concentrations and 
for k near the crossing of the dispersion with the shadow 
of the FS near (7r/2,7r/2) (a), for k at the FS near (7r,0) 
(b) and for k at the shadow of the FS near (tt, 0) (c- 
d). Qualitatively, these data agree with our findings for 
the previously used model compound since we find pro- 
nounced occupied and unoccupied shadow states in the 
BZ as small additional satellites in the spectral density. 
Again, these states are not related to new quasi particles 
but rather to an incoherent amount of spectral weight 
caused by an increased scattering rate [Q. However by 
comparing the data for x = 0.02 and x — 0.08 as in Fig. 
2 (c) and (d), we see that they become much weaker due 
less pronounced nesting of the FS and the corresponding 
much smaller magnetic interaction Vq{uj) as can be seen 
in the inset and in Fig. 1. Comparing these results with 
the t' =t" =0 FS we find that for x = 0.08 the shadows 
are rather strong Thus the observation by Aebi et 
at ^ can not be satisfactorily understood by considering 
only a single Cu02 plane which is in agreement with the 
study of Ref. |^ . since the shadows were found near the 
optimal doping at x ~ 0.15. 

However, as recently discussed for the t' = t" = 
model dispersion [E0| the consideration of a finite t± in- 



creases the shadow state intensity due to the additional 
bilayer correlations. Here the fact that for low doping 
VqL(w) w V^{uj) lead via Eq. | to V++{uj) w and 

(lu) « Vq (uj) such that we find that the spectral 
weight is not only shifted by the momentum Q, but 
simultaneously also from the bonding to the antibond- 
ing band and vise versa. Now by taking a finite t±_ 
into account we find for YBCO and BSCCO-like sys- 
tems that the inter-layer antiferromagnetic coupling also 
increases the shadow state intensity and they start to ap- 
pear for X < 0.12 with a maximum intensity at tj_ ~ OAt. 
Furthermore as demonstrated in Fig. 3 for YBCO we 
find that the most favorable region to observe shadow 
states in ARPES is in the neighborhood of the {tt/2, tt/2) 
point, where main and shadow band are well separated. 
Near (tt, 0) the absolute intensity of the shadow states 
is largest, but they are difficult to detect because of the 
superposition of shadow peaks and the dominant main 
band contributions. Furthermore the strongest shadow 
band signal might not be obtained with photoemission 
but with an inverse photoemission measurements (IPES) 
as shown in Fig. 3 (a) where a pronounced shadow 
peak appears. However the current limited resolution 
of IPES measurements might impede the observability 
of these states. For the BSCCO system where we use 
a k-dependent inter-layer hopping tj_ (k) as suggested be 
band structure calculations p^ , we find similar results 
near (tt, 0). However on the diagonal at {tt/2, tt/2) both 
bands are not splitted and the shadows of the bonding 
and antibonding band are at the same position. 

The overall shape of the quasi particle dispersion in the 
bilayer compound YBCO and the shadow band formation 
is presented in Fig. 4 where we focus our attention on the 
bilayer splitting. The experimental findings concerning 
the bilayer splitting in the cuprates are still controversial. 
For YBCO Liu et al. [|] and Gofron et al. [|j found indi- 
cations via ARPES for a spfitting A^^P£(7r, 0) = 110 meV 
that is much smaller that the theoretical predictions from 
band structure calculations jl^] whereby for BSCCO dif- 
ferent photoemission groups come to different conclusions 
concerning the existence or the absence of a finite bilayer 
splitting I^J^. However, from our results for the bilayer 
Hubbard model we conclude that the strong antiferro- 
magnetic correlations in the high-Tc superconductors are 
responsible for the smallness of the bilayer splitting and 
might explain these interesting experimental results as 
also discussed by Liechtenstein et al. |2^ . 

At this stage, it is of interest to compare our results 
with the interesting argumentation of Vilk that 
shadow states and a variety of related phenomena ap- 
pear only if the magnetic correlation length ^ is larger 
than the thermal dc Broglie wavelength Af^ = vp /ttUbT . 
Here vp \s the Fermi velocity of the system. Assuming 
that vp is given by its uncorrelated value vp ~ ttu t (a: 
lattice constant) would yield shadow states for low tem- 
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peratures only for extremely large magnetic correlation 
length. However, our self consistent calculations yields 
a self stabilization of the shadow band phenomena: the 
increasing quasi particle scattering rate which is related 
to the formation of shadow states gives rise to a pro- 
nounced flattening of the correlated dispersion as can be 
seen in Fig. 4. This causes a renormalization and a de- 
crease of the Fermi velocity for decreasing T. This mech- 
anism allows observable shadow bands even for low tem- 
peratures since Xth remains almost constant as a function 
of T and consequently Xth < £. can be fulfilled for moder- 
ate correlation length as observed experimentally in the 
high- Tc superconductors. Note that the limit T ^ is 
still an open question although we believe that for very 
low temperatures the susceptibility and the correlation 
length ^ behave differently compared to the effective in- 
teraction in the expression for the self energy ||2^. Thus 
this demonstrates impressively the importance of a self 
consistent calculation of the Fermi velocity that decreases 
down to T « 40 K by taking explicitly the changes of the 
spin fluctuation and the quasi particle spectrum into ac- 
count. 

In Fig. 5 we present our data for the superconduct- 
ing order parameter for the YBCO model where we plot 
the frequency dependence of the superconducting order 
parameter (/)k(ijj) at k = (tt, 0) for x — 0.08. Here we 
treated the mono and bilayer Hubbard model below Tc by 
solving the strong coupling Eliashberg equations thereby 
considering the full momentum, frequency and doping de- 
pendence of the intra- and inter-layer interactions. Here 
one obtains one order parameter for each band, namely 
(j>^{uj). These are connected to the layer representation 

via = ± cjj^iuj), where ^(w) de- 

scribes intra- (inter-) layer Cooper pair formation. By 
solving these equations for the YBCO dispersion with 
tj_ = 0, we find for all x a d.j.2^y2 wave superconducting 
state (Tc = 70 K for x — 0.08) without the need to in- 
troduce phenomenological interactions that enforce the 
formation of Cooper pairs |24j2^. Note that this value 
is only slightly smaller than the corresponding value for 
t' = t" = {Tc = 97 K for x = 0.12) whose FS is rather 
similar to the La2-3;Sr3;Cu04 (LSCO) compound. Con- 
cerning the pairing symmetry in a bilayer system also find 
dx2_y2 pairing symmetry and similar values for the tran- 
sition temperatures. Interestingly, as can be observed in 
Fig. 5 the d-wave state is characterized by an increasing 
contribution of coherent inter-layer pairing with decreas- 
ing doping, where for x — 0.08 Cooper pairs are formed 
by electrons from the same layer as from different layers 
with almost equal probability. 

Although it is a remarkable success that spin fluctu- 
ation induced pairing interaction cause superconductiv- 
ity in mono and bilayer compounds, it might be even 
more important to explain the rather different transition 
temperatures in the large family of the high- Tc super- 



conductors. Here a comparison of our monolayer results 
for LSCO with the data of the bilayer YBCO compound 
shows that Tc of YBCO should be even slightly lower 
than for LSCO the system when all other parameters like 
U or the hopping t are taken to be the same which is def- 
initely in contradiction to the experimental observations 
(yysco ^ 90 K versus Tc^^^-o 30 K). 

However it has been recently argued that the larger Tc 
is YBCO is mainly due to finite size effect and to the rel- 
atively strong coupling of the unit cells ||2^ , which is not 
included in our approach. This point is also supported 
by experiments on ultrathin YBCO layers separated by 
insulating PrBa2Cu307_5 which show that one separated 
bilayer of YBCO yields Tc « 20 K and that the coupling 
of these bilayers boost the transition temperature to its 
bulk value p^ , ^^ . In this context our findings that a 
single bilayer of YBCO has a smaller Tc that LSCO are 
in qualitative agreement with the experiment although 
there is still the need for a consideration of inter-bilayer 
effects and for understanding the importance of other in- 
gredients of the high- Tc materials for superconductivity 

In conclusion, we investigated the doping dependence 
of the mono- and bilayer Hubbard for the YBCO and 
BSCCO systems by using the FLEX approximation. 
The strong antiferromagnetic coupling between the lay- 
ers leads to an enhanced shadow state intensity and to 
their observability in YBCO- and BSCCO-like models. 
We found that the bilayer splitting is reduced 50%) 
and the inter-layer hopping is effectively blocked for small 
doping. Finally, we presented results for the supercon- 
ducting which has a d-wave order parameter for mono- 
and bilayer compounds whereby the later one has signif- 
icant inter-layer Cooper pair contributions. 
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FIG. 1. Effective interaction Vc^{^^ ~ 0) of the BSCCO/YBCO model. Note the commensurable structures, e.g. the fact that 
Vq(ijJ — 0) is peaked at (tt, tt) as observed in experiments. 
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FIG. 2. Spectral density pkicj) for the BSCCO/YBCO dispersion for k points as indicated in the inset and two doping 
concentrations, (a): Crossing of the shadow of the Fermi surface on the diagonal in the BZ. (b): Fermi surface crossing near 
(tt, 0). (c) and (d): Crossing of the shadow of the Fermi surface near (tt, 0) for x = 0.02 and x = 0.08. 
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FIG. 3. Formation of shadow states for YBCO. (a): p^{to) near (7r,0). (b): p^{iij) on the diagonal near (7r/2,7r/2). For 
simplicity only the FS of the bonding band and the shadow of the antibonding band are shown in the inset. 
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FIG. 4. Quasi particle dispersion of the bilayer BSCCO and the YBCO model (solid lines: main bands, dashed lines: shadow 
bands). Note the strong suppression of the bilayer splitting from Ae = 200 meV for (7 = to Aej.^ g) = 86 meV. 
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FIG. 5. Superconducting order parameter Re (/)k(<j-') for the monolayer {t± = 0) and for the bilayer with 
k — (n, 0) and the YBCO dispersion. Note the large contribution of the inter-layer Copper pair formation. 
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